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Abstract 

We solve the B.S. equation of glueballs under instantaneous ap- 
proximation. With the B.S wave function obtained we calculate 
the decay width of glueballs to two pseudoscalar mesons, r(7T7r), 
T(kk) and r(i]f]). Y{rfq') from QCD anomaly is also estimated. 
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1. Introduction The existence of the glueball states is the prediction of QCD, whose 
discovery would be a direct confirmation of the non-ablian character of QCD. A lot 
of investigation both in theory [1-5] and in experiment [6-8] have been made in the 
glueball searches. At present there are several experimental glueball candidates, such 
as / (1500), /j(1700) and £(2230), but non is identified as glueball definitely. It's 
difficult to identify a glueball because so far we haven't understood the glueball clearly 
in theory: we don't know clearly its mass, its production rates and its decay width. 
Therefore, more investigations are necessary. 

In this paper, we focus on ++ glueball's decays. Up to now, / (1500) is considered 
as the lowest-lying ++ glueball candidate. One of the most important reason is that 
the branching ratio of decay /o(1500) — > rjrj' is surprisingly large, i.e T(j]T]')/V(t]T]') ~ 
1. Since not only 7777' final state is almost on the kinematic threshold but also it is 
suppressed by SU(3) flavor symmetry, the naive quark model gives r(rp]')/T(rp]') ~ 0. 
Therefore if must be produced by another mechanism. One of the possible mechanism 
is the famous QCD anomaly. However, it is well known that QCD anomaly is large N c 
suppressed, in order to produce enough rate of 77' production, / (1500) should have a 
high gluon density, i.e., it is a glueball. On this assumption , the further theoretical 
calculation is still necessary to check the rate r(r]i]')/r(r]rj) and T(rjr]')/r(7CTT). We 
notice there have already been some works on this topic from lattice QCD approach, 
however, an model-dependent calculation is still useful and may give some helpful hints. 

In an early paper ||, we studied the glueball spectrum in the framework of B.S. 
equation, and the results are in good agreement with that of Lattice calculation. This 
encourages us to do more investigation. By solving the B.S. equation about gluon- 
gluon bound state, we can obtain not only the glueball spectrum but also the glueball 
B.S. wave function. In this paper we try to use these wave functions to estimate the 
decay width of glueballs to light pseudoscalar mesons in the framework of perturbative 
QCD. 

This paper is arranged as follows. In the next section we briefly review the construc- 
tion of the B.S. equation for gluons bound state. In section three, first we calculate 
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the glueball width to two pions. Then we give the estimations of width to two K's and 
two 77's. In section four, we give the estimation of T(f]i]') from QCD anomaly. In the 
last section a discussion is presented. 

2. The Glueball B.S. Equation The details of the construction of the glueball B.S. 
equation has been given in ref . || . Here we only present a brief review. 

Let A a (x\) and A b (x 2 ) be the gluon fields at points x\ and x 2 , \G) the bound state 
of two gluons with mass Mq and momentum P M . Then the B.S. wave function for a 
bound state is defined as 

x£(P,*i,a: a ) = (0\T(A«( Xl )A b u (x 2 ))\G), (1) 

where //, v are Lorentz indices, and a, b color indices. The glueballs are color singlet 

states, so we have 

X %(P,x 1 ,x 2 )=6 ab X ^(P,x 1 ,x 2 ). (2) 
With a standard method we obtain the B.S. equation for a color singlet glueball state, 

X»„(P, q) = A m ( Pl )A^(p 2 ) J 0^G a ^(P, q, k) Xf „(P, k), (3) 

where X^(PiQ) is the Fourier transformation of X ^ U (P, x\,x 2 ). 

The tensor kernel G af3pr7 (P, q, k), in which the color indices are suppressed, is defined 
as the sum of all two-particle irreducible graphs, and A^pi) is the full propagator of 
the constituent gluons with momentum p^. 

P = Pi+P2, 2q = p 1 -p 2 , (4) 

Under instantaneous approximation, we ignore the go an d ko dependence of the 
kernel G al3pa (P,q,k). Further we replace A^pi) with a free propagator. Then in 
Coulomb gauge we obtain a three dimensional equation 

E(M 2 G -AE 2 )^(P, q ) = 

i{5u, - q -^)(S jf - ^) / 70 G ^(^ q, k)MP, k). (5) 
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where E = y/q 2 + m 2 {m is the constituent mass of gluons, which will be chosen as 
zero in the present calculations), and <Pij(i,j = 1,2,3) is the three dimensional wave 
function, which is related to Xij by the following equation 

W(P,k) = J dk oX ^(P,k). (6) 

The kernel is divided into two parts: the short distance part G^ and the long dis- 
tance part Gfj kl . G\j\ t is obtained approximately by calculating the three lowest order 
diagrams shown in fig.l. 
Besides equations (4), we also have 

P = P3+P4, 2k = p 3 -p 4 . (7) 

Calculating the diagrams a, b and c shown in fig.l, Gfy pa (P, q, k) can be expressed 
explicitly as 

" rO r'O ri r'j j , 

+ * L fi-O y «-ir) 



l 2 I 2 v 13 \ 2 



- (2g^ v g pa - g w g V(J - g^9u P ) J , (8) 
where /(= q — k) is the momentum exchanged between the two constituent gluons, and 

CppAPuPs) = (Pi - 2p3)p9 P T + (Pi + P3)r9p P + (ps - 2pi)p9pT, (9) 

C uaT/ (p 2 ,P4) = (P2 - 2p4)vgaT> + (P2 + Pi)r'9ua + {Pi ~ 1p<l)a9vT'- (10) 

The factor 3 on the right-hand side of equation (8) is the color factor which is | in 
the case of qq bound state. Diagram a and diagram b make the same contribution to 
physical states, so there is a factor 2 in equation (8). The strong coupling constant a s 
is chosen as a running one, 

127T 1 

®s = ^— — ; 13 — r, (11) 

27 l n (a + W — ) 

IV QCD 

where a is a parameter introduced to avoid the infrared divergence. We take a=4.0, 
which implies that the largest value of the running coupling constant is 1.0, and we 
choose Kqcd = 200MeV. 



As for the long distance part Gfj kl , we only choose in the three-gluon-vertex (9,10) 
the terms containing tensor g w g va as the spin dependence of the confining part, because 
such terms have nothing to do with spin effect. Therefore we assume 

G { X a = 2i{pi + Ps) ■ {P2 + P4)g w 9uaG{l), (12) 

where G(l) is spatial dependence of confining part of the kernel. We choose 

G(f) = ?£, (13) 

which corresponds to a linearly growing potential. A is the string tension, which 
in quark potential model takes the value 0.18(GeV) 2 . However, here we chose A = 
0.36(GeV^) 2 . We expect that confining force between two gluons is one times stronger 
than that between two quarks because each gluon acts like a qq pair [4]. Of course, 
such a potential is very singular at the point I = 0. Some form of regularization is 
necessary. Our method is given in ref.f|. For ++ glueball, the B.S. wave function can 
be expressed as 

^(P,q) = / S (q)(^-^). (14) 

With the kernel and the form of the wave functions chosen as above, the equation can 
be solved numerically. The glueball spectrum obtained are in good agreement with 
that of lattice calculation. At the same time the wave functions (14) and (15) are 
also obtained. Such wave functions should be normalized before being used in the 
calculation of the decays. If we neglect the P$ dependence of the kernel, then the three 
dimensional wave function fij(P, k) is subjected to the normalization condition 

/ 7?^Mq)w«(q)(*tt - q -^){h - q -§) = (2vr) 2 (15) 

3. The Width of ++ Glueball to Two Liglt Mesons 

i. The decay to two pions In the framework of perturbative QCD, the decays 
of the glueball to two pions can be expressed by the diagram shown in fig.2. The 
amplitude is expressed as 

M(G(0 ++ ) - 7T+7T-) = J^d^d&tpA&M&iMuky + fci • M/.h 

(16) 



where a s is effective coupling constant which will be discussed below, k\, k 2 are the 
momentum of pions. is the pion wave function or distribution amplitude. Xij 

is the four dimensional glueball B.S. wave function defined in section 1., which can be 
replaced by the three dimensional wave function cp^ through 

|q|(MS-4q») (17) 

(fij is obtained by solving the B.S. equation. 

The mass of pions can be neglected, and thus we have 



h ■ k 2 



G 
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ii 16 + 61 , , / in >. 

|q| = M G . (19) 

With these relations, equation (16) is reduced to 

Aja Mr<f 2 r 1 

M (G(0 ++ ) - 7T + 7T-) = S 9 GJ * ^6<^(6K(6) 

There are several forms of the pion wave function. Here we will use the form given by 
Chernyak by QCD sum rule||. 

Mt)=^?0--?)- ( 21 ) 

In equations (16) and (20), a s is the effective coupling constant. In fact, we should 
have taken ^Ja s ■ a s at the points corresponding to the gloun virtualities: ct s {p\) = 
a s ( — — — — -—^M G ) and a s {p\) = a s ( — — — ~^ 2 M G ) . However, for simplicity, we 
replace pi 2 by their character values pi 2 , and then extract a s out of the integrals over 
6 and 6- Pi 2 can be determined by the forms of pion wave functions. Therefore, d s in 
formula (16) is chosen as 0ls ^ 1 }^2^\ • ^ or Chernyak's form, the pion wave function 



a s {M G /A) 

iximum at £=0.7, ; 

glueball 



have maximum at £=0.7, and p\^ = — — — 2 M G . Therefore, a s = 0.69 for 0^ 
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Substituting expression (14) into formula (20), we get 

8ia,fiM G I 



M(G(0 ++ ) -> 7T + 7T-) = 



9 



(22) 



with 



= y_ 1 ^1^2^(6)^(6) — n _£2)n _ — Ml**!) 



(i-a(i-a 

Finally, after making the phase space integration, we obtain 

4a 2 s M G f%P 



r(G(0 ++ ) -> 7T + 7T-) 



817T 



(23) 



The numerical results is shown in table 1. 

ii. The decays to two K 's and two r\ 's To estimate the glueball decays to two k 
and two 77 mesons, we construct the following chiral Lagrangian in the lowest order 



L eff = gct>tr{d^ + d^) 



where is the field of the scalar meson, and £ = e i* , 



(24) 



/ ^0 



7T = 



V 



7T" Tj 

k~ 



7T° 77 

72 76 



2r] 



k° --7= 



>/6 / 



g is a coupling constant which can be determined from the result of two pions decay 
obtained above. Neglecting the factor of the phase space, we have 

r(7r + 7r-) = r(A; + Ar) = r(k°k°) = ^(ttV) = ^r{ vv ) 

Taking the factor of the phase space into consideration, We give the numerical results 
in table. 1. 

4. The decay of glueball to r\, rj' form QCD anomaly As mentioned in the introduc- 
tion, 77, v( final state cannot be produced enough only from naive quark model, there 
should be another mechanism to produce them, namely the QCD anomaly (see fig. 3) 

3 a. 



An 



m q q^q 



(25) 



q=uds 



where j£ is the singlet axial current [10]. Indeed it is the anomalous contribution to 
the divergence of the axial current that makes the 577(3) singlet state much heavier 
than the Octet pseudo-Goldstone state and distinguish 77 from rj'. Obviously, keeping 
SU(3)l <E> SU(3)r flavor symmetry unbroken, 77 cannot be produced from two gluons. 
Actually, physical 77, rf states are the mixing states of flavor eigenstate rji and r)$ 

77 = cos9 v r] 8 — sin8 v T]i, (26) 

7/ = sin9 v r]s + cos9 v r)\, (27) 

where the mixing angle 8 V is estimated in the range of —10° to —20°, so rj' is dominantly 
771. The effective coupling of r\\gg can be written in the form 

H(ql q 2 2 , ql^e^qU^/cosO, (28) 

where qx, q2 and e\, e 2 are the 4-momenta and polarizations of two gluons and a, b 
are color indices. H is a form factor that is a general function of momenta, qf, q\ 
and g 2 ,. For simplicity, we set H as a constant , then, from the experimental data 
T(J/Psi -> 77/)/r(J/tf -> e + e _ ) = 7 x 10" 2 , one may obtain H ~ 1.8 GeV" 1 Q. 
Corresponding to fig. 3, the amplitude dominated by on-shell gluons can be written as 

d 4 q cPqx d 3 q 2 d 3 q 3 d 3 q 4 



M(G -> 7]7]') = 3V8cos6 v sin6 v J 



(2tt) 4 (27r) 3 2| Pl | (2tt) 3 2| P2 | (27r) 3 2|p 3 | (27r) 3 2|p 4 | 



■(27T) 4 5 4 (ip + q- qi - q 2 )(2n) 4 5\k 1 q z ){2ix)H\k 2 - q 2 - q A ) 

■ X ,u(P, q)C^(- gi , q 2 )C a ^(-q 3 , q^e^e^q'q^ql (29) 

where C^(-gi, q 2 ) and C a/3u (-q 3 , q 4 ) are the three-gluon- vertex defined as eq.(9), and 
3^ is a color factor. With the the glueball B.S. wave function obtained in section 2, 
the calculation of the amplitude is straight forward. The corresponding width is list in 
table 1. 

Because the sum of and m v i is quite near the mass of glueball, M, we expect the 
above result from QCD anomaly will depend on M. On the other hand, the glueball 
mass is not exactly known either from experiment or from theory. Therefore, we show 



the dependence of Ta^ m i on glueball mass M in fig. 4. The value of T(rjr},) in table 1 
corresponds to M — 1.6GeV. 

5. Discussion and Conclusions 
In this paper, we study the decays of ++ glueball to two light mesons by using the 
constituent gluon model. The glueball wave function is obtained by solving B-S equa- 
tion, the pion is light so that it can described by the light cone wave function. SU(3) 
chiral symmetry is used to predict decay widths T(KK) and F(j]T]). we also study the 
i] rf production from QCD anomaly, which almost is zero from the naive quark model 
prediction. We are more interested in the ratios of the different decay channel, since 
the ratios is less sensitive to the glueball wave function we use. From table 1, we get 

T(t]T]')/T(t]T]) = 0.82 r(vr 7r )/r( W ) = 1.5 



Such ratios have values in experiment for / (1500). The Crystal Barrel Collab. [|i~2]| 
gives 

r(r)r}')/r(r)r)) = 0.29 ± 0.1 (30) 

r(7r°7T ) /T(r)T]) = 1.45 ± 0.61(3tt° channel) (31) 

r(7r°7r°)/r( W ) = 2.12 ± 0.81 (Coupled - channel), (32) 



GAM group gave a rather large ratio jl3" | 



r (W)/ r (w) = 2.7 ±0.8 (33) 
We see that our results agree with the experiment data of /o(1500). 
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Figuer Caption 

Fig.l: The diagrams contributing to the short distance part of the B.S. kernel. 
Fig. 2: The decays of glueball to two pions in purterbative framework. 
Fig. 3: The decays of glueball to r\r( from QCD anomaly 
Fig. 4: The dependence of T(rjr}') on glueball mass. 
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Table Caption 



Table 1: The partial width of the glueball to two light mesons (MeV). 
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Table 1. The partial width of the glueballs to two and four pions(MeV). 
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